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Low energy elastic scattering between clouds of Bose condensed atoms leads to the well known 
s-wave halo with atoms emerging in all directions from the collision zone. In this paper we discuss 
the emergence of Hanbury Brown and Twiss coincidences between atoms scattered in nearly parallel 
directions. We develop a simple model that explains the observations in terms of an interference 
involving two pairs of atoms each associated with the elementary s wave scattering process. 

PACS numbers: 03.75 Gg, 34.50-s 



I. INTRODUCTION 

In a number of experiments, Bose-Einstein condensates 
have been prepared to collide with each other with well 
defined collision energies and momenta. At the micro- 
scopic level, when two particles with equal mass and op- 
posite velocities collide in an s-wave collision, the colli- 
sion partners will propagate away from each other with 
the same probability amplitude in all directions, but 
their individual momenta are correlated in opposite di- 
rections, as their total center-of-mass momentum is con- 
served in the collision process. In experiments with col- 
liding condensates, the scattering into all directions has 
been clearly observed as a so-called s-wave halo of scat- 
tered particles The observation of pair correlations 
of particle leaving the collision region back-to-back, see 
Fig. 1, requires efficient detection of all momentum com- 
ponents of individual atoms, and this correlation has re- 
cently been observed as a significant coincidence signal in 
a collision experiment with Bose-Einstein condensates of 
metastable atomic helium [2] ■ The same experiment also 
observed an increased coincidence of particles scattered 
into the same direction. This phenomenon is due to the 
bosonic nature of the particles and to the fact that several 
independent scattering processes occur simultaneously. 

We shall present a theoretical analysis of this Hanbury 
Brown and Twiss correlation phenomenon, aiming at a 
simple model which explains its qualitative and quanti- 
tative character in experiments. It is important to em- 
phasize that the appearance of atoms moving in the same 
direction is not compatible with momentum conservation 
in a single collision process of two counterpropagating 
atoms, and our discussion will, indeed, refer to effects 
that rely on a many-body treatment of the collision of 
larger ensembles. In order to get physical insight, we 
will separate the problem in two: first, we treat indepen- 
dent pairwise collisions in which the coUisional interac- 
tion gives rise to pairs of scattered atoms, and second, we 
treat the evolution of the many-body state describing the 
ensemble of scattered atoms neglecting the interactions 




FIG. 1: Diagrammatic representation of two condensates col- 
liding and giving rise to an s-wave halo of scattered particles. 
Particles are scattered pairwise back-to-back. 



in order to apply analytic methods. The validity of this 
separation and means to improve the theory, if necessary, 
will be discussed. 

In Sec. II, we present a full second quantized descrip- 
tion of the collision of a large number of identical bosons. 
We shall write down the second quantized many-body 
Hamiltonian, and discuss how the elementary processes 
of interest relate to the different terms in this Hamil- 
tonian. We treat the case where all the atoms initially 
populate two counter propagating single particle states 
which are only weakly depleted by the collisions, and 
which will hence serve as c-number field sources for cre- 
ation of pairs This is in analogy with the quantum optics 
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treatment of spontaneous four-wave mixing, where pho- 
ton pairs are generated from the interaction of two inci- 
dent laser beams, described by classical electromagnetic 
waves. In Sec. Ill, we discuss the Bogoliubov transfor- 
mation which provides a very accurate approximation to 
the time evolution of the system. We shall not, how- 
ever, apply this transformation in a quantitative treat- 
ment, but rather show that its formal structure already 
predicts the coUinear (Hanbury Brown and Twiss) cor- 
relation and motivates a quite general analytical Ansatz 
for the quantum state of the scattered atoms. In Sec. IV, 
we shall consider the leading two- and four-atom terms 
in an expansion of the quantum state of scattered atoms, 
and show that they hold the key to the observed Han- 
bury Brown and Twiss correlations. In Sec. V, we shall 
use energy conservation and phase matching considera- 
tions to motivate a simple analytical model, from which 
we show that the coincidence of scattered particles in the 
same direction, though a many-body effect, can be un- 
derstood quantitatively from the properties of the simple 
two-atom scattering wave function. Sec. VI concludes 
the paper with a discussion of the insights offered by our 
analysis. 



II. COLLIDING BOSE-EINSTEIN 
CONDENSATES 

The Hamiltonian, 

H = J dfii^f) (^-^A + V{r) + |^t(f)*(f)^ #(f), 

(1) 

with field operators obeying the bosonic commutator re- 
lations [^'(r), ^'^^(r')] = 5{f — r'), gives a good descrip- 
tion of bosons interacting at low collision energies via a 
short range potential. The interaction is represented by a 
delta-function interaction term with strength g, propor- 
tional to the s-wave scattering length. The atoms may 
be subject to a wide range of trapping or guiding poten- 
tials V{f), or they may propagate freely (V = 0), and 
the initial state of the system may be specified according 
to experimental preparation procedures to describe for 
example a single condensate or several macroscopically 
populated components. We are interested in the situa- 
tion, where two condensates with well defined opposite 
momenta, and hence with relatively large spatial extent 
propagate towards each other. The conventional second 
quantized Hamiltonian fully describes the problem, and 
a Monte-Carlo type simulation of the dynamics [1, |j, Q j 
and perhaps even simpler simulation approaches based on 
truncated Wigner function expansions [gI, I?! , may solve 
this problem in full generality by full 3D propagation of 
stochastic Schrodinger type equations. 

We assume that elastic collisions occur with a suffi- 
ciently small cross section that the colliding condensates 
are only weakly depleted due to the collision term in ([1]) . 
The Hamiltonian has terms describing the kinetic energy 



and the potential energy of atoms moving in the external 
potential and finally, a term describing the mean field 
repulsive or attractive potential due to the other atoms 
of the colliding condensates. But the product of two cre- 
ation and two annihilation operators in the interaction 
term does not only read as density dependent correction 
to the potential energy in the Gross-Pitaevskii equation: 
the product of two creation operators may also cause the 
creation of a pair of atoms with momenta entirely differ- 
ent from the incident ones, extracted consistently from 
the condensates by the product of annihilation operators. 
The pairs of atoms "created" in the scattering process are 
the ones that are detected as the s-wave halo around the 
condensate collision region in Fig. 1. 

We can think of each point in the collision zone as a 
point source for a pair of initially close atoms (atoms 
only collide at short range), which are subsequently sep- 
arated by free propagation, perturbed by the interac- 
tion with the condensate components. This propagation, 
together with the coherent addition of pair amplitudes 
originating from the entire collision zone leads to a com- 
plicated many-body entangled state, but energy conser- 
vation, imposed after sufficiently long interaction time, 
and momentum conservation, imposed by phase match- 
ing, serves to justify our simpler model, described below. 



III. BOGOLIUBOV APPROXIMATION 

If the original condensates are only weakly depleted 
by the scattering, we may expand the Hamiltonian to 
second order in operator terms and get linear Heisenberg 
equations of motion that may be solved and tell precisely 
which Bogoliubov transformation exactly solves the prob- 
lem. This approach was followed in Refs. The 
Bogoliubov solution yields an expression for the atomic 
annihilation operators at any given time, expressed as a 
linear combination of the annihilation and creation oper- 
ators at time zero, where the initial state is assumed to 
be known (incident condensate wave functions, no scat- 
tered atoms). The mean atom number and any higher 
order correlation function of the field can therefore be 
expressed algebraically in terms of the expansion coeffi- 
cients of the Bogoliubov transformation and the known 
vacuum expectation values of field operator products. 

Even though the full many problem has thus been 
reduced to partial differential equations of a complex- 
ity comparable to the single particle Schrodinger equa- 
tion, one still has to solve time dependent wave equations 
in three spatial dimensions. Here we shall demonstrate 
some properties of the solution that follow by purely an- 
alytical arguments, i.e., without access to the precise so- 
lution. 

Although obviously related, the use of the Bogoliubov 
transformation here is different from the Bogoliubov ap- 
proximation used to identify low-lying, collective excita- 
tion modes in a condensate. The analysis rather follows 
the philosophy of squeezed light generation with opti- 
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cal parametric oscillators in quantum optics, where the 
Bogoliubov method is used to diagonalize a multi-mode 
Hamiltonian with pair creation and annihilation opera- 
tors [Ui], see also [ll, HI]. 

Since the initial state of the atomic scattering modes 
(momentum components) of interest is the vacuum state, 
which has a Gaussian (Wigner) probability distribution 
for the multi-mode field variables, and the Bogoliubov 
transformation is linear in field operators, the state will, 
independently of the precise form of the transformation, 
at all later times be a Gaussian with vanishing mean field 
expectation value [l^ . If we restrict the analysis to a sin- 
gle final momentum state (mode), by a partial trace over 
all other modes, the state of this mode is also a Gaussian 
state with vanishing mean amplitude. It is thus fully 
characterized by the second moments of the hermitian 

linear combinations = (^^\k) + '^(k)^ /\/2, pj: = 

i (^^'t(fc) - *(A:)^ /V2 of the field operators. We now 

wish to establish that our Gaussian distribution is sym- 
metric, i.e., Var(q^)=Var(pj:). This indeed follows if the 
"anomalous" moments (^t(fc)2^ ^ (*(A:)^) = 0, i.e., if 
there is no coherence between states differing by two 
atoms propagating in the given direction. We now ap- 
ply the physical argument, that the coUisional Hamilto- 
nian does not produce such coherence, since the collision 
process can only produce pairs of atoms propagating in 
opposite directions, and states, e.g, with zero and two 
atoms with momentum k must also contain zero and two 
atoms with momentum around —k. The anomalous mo- 
ments vanish due to the orthogonality of these parts of 
the wave function. 

It is well known in quantum optics, that a symmetric 
Gaussian state is equivalent to an incoherent mixture of 
number states with exponential number distribution, also 
known as a thermal state with the density matrix [11], 

Pi = {l-mJ2\tr\n){n\. (2) 

The state conditioned upon detection, and annihilation, 
of a single particle reads, 

where is a normalization constant. A straightforward 
calculation shows that this state has precisely twice as 
many bosons on average as ([2]) , and hence that the prob- 
ability to detect two bosons by a low efficiency detec- 
tor is twice the square of the single quantum detec- 
tion probability. It thus follows that the coincidence 
counting rate for observing two atoms leaving the col- 
lision zone in the same, narrowly defined, direction, 
(#t(fc)\I't(fc)\p(fc)\i/(fc)) is twice the square of the mean 
counting rate, and twice the coincidence rate for seeing 
atoms in two unrelated directions. 

Without performing any calculations, we therefore 
understand qualitatively the observed coincidences ob- 
served in the experiments Q as the direct consequence 



of the thermal counting statistics (Gaussian quadrature 
distribution) of the output flux in all scattering direc- 
tions. This is the famous Hanbury Brown and Twiss ef- 
fect [13, EM 0, [ioj observed originally as photon bunch- 
ing in chaotic light resulting from the addition of the 
contributions of many incoherent emitters. In order to 
provide a natural estimate of the HBT momentum cor- 
relation function, one could develop the field correlations 
by solution of the linear Bogoliubov-de-Gennes equations 
for the problem [1, [^, [13] which by the corresponding 
linear transformation of operators provides the first and 
second order momentum correlation functions and hence 
the momentum range within which the bunching effect 
takes place. Here, we will rather keep track of the binary 
scattering states, and in particular of the counterprop- 
agating partners, which will give us an alternative and 
very useful physical interpretation of the effect. 

The Bogoliubov transformation of field operators is 
equivalent to a multi-mode unitary squeezing operation 
[111 ], which is indeed nothing but the time evolution op- 
erator of a Hamiltonian with quadratic terms in field cre- 
ation and annihilation operators. Such an operator can 
be ordered as a product of three exponentials [13, [l3| : 
one involving a sum of products of pairs of creation op- 
erators, one involving a sum of products of creation and 
annihilation operators and one involving a sum of prod- 
ucts of pairs of annihilation operators. When acting on 
the initial vacuum state vector, only the unit term of 
the series expansion of the latter two exponentials con- 
tribute, and the state can therefore be written in terms 
of a quadratic form of creation operators of atoms, e.g., 
in the momentum space representation, 

]* A^*exp( j dkidk2^/j{ki,k2)'l>\ki)¥ {k2))\vac> . 

The function ^(A:i,/c2) generally depends in a complex 
manner on the dynamical evolution. It is of course re- 
lated to the scattering wave function of a single pair of 
atoms, and we shall come back to this relationship in con- 
nection with the model studied in Sec. V of the paper. 
The use of second quantization automatically yields the 
bosonic symmetry of our state, but in addition we can 
require that the pair amplitude function obey the explicit 
exchange symmetry 'ip(ki,k2) — '4'{k2,ki). For now, let 
us assume, that the propensity for atoms to be scattered 
into opposite direction also implies that V'(^i,fe) takes 
non-vanishing values for all directions of the scattered 
particles, but only if ki ~ — fc2- The fmiction ijj(ki,k2) 
is not a normalized wave function: the larger its ampli- 
tude the more particle pairs are created, and higher order 
terms of the exponential play more and more important 
roles. The many-body state \'^) is normalized by the 
prefactor in 

We now proceed to determine the density-density cor- 
relations of atoms detected in two different directions, 
labeled by momentum states (fc, fc'), i.e., the expectation 
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value 

F{k, k') cx ('I'+(fc)*+(fc')*(fc')*(fc)> (5) 

The state ([4]) is a Gaussian state, and by Wick's the- 
orem [2l| this expectation value can be written down in 
terms of only pair-expectation values. We shall address 
the contribution from the four-atom component in the 
expansion of the exponential in (j4|), as this provides a 
straightforward interpretation of the origin and the be- 
havior of the atomic Hanbury Brown and Twiss correla- 
tions. 

IV. TWO-ATOM AND FOUR-ATOM STATES 

The state ^ can be written explicitly, 

|*> ^ iV^,(|mc> 

-I- J dkidk2^p{ki,k2)i''^{ki)ii+{k2)\vac > 

+ ^{j dkidk2i^{kiM)'^^{ki)-^^{k2)f\vac> 

+ ...). (6) 

The zero order term is the vacuum state. The first order 
term is a two-atom state of atoms propagating back-to- 
back, and the second order term of the series expansion 
of ([4]) is the four-atom state 

1*4) = ( I dk,dk2tP{k,,k2)i'+{ka)^+(k2))^\vac >, (7) 

which we will show accounts for the observed HBT ef- 
fect. The squared pair creation operator in ([7]) can be 
expanded as a four-fold integral. To obtain the correla- 
tion function we have to apply the product of the two 
annihilation operators on |\E'4) and determine the squared 
norm of the resulting state, 

F{k,k') (x\\i>{k)i>{k')\^4W. (8) 

Using the field commutator relations, we can shift the an- 
nihilation operators to the right of all creation operators 
in ([H]). This yields a total of 12 terms, which by relabel- 
ing and use of the exchange symmetry can be reduced to 
a sum of three different contributions, 

*(fc)*(fc')l*4) OC 
J dhdk2Mki,k2)yj{k,k') + i:{ki,k)^/j{k2,k') 

+%lj{ki,k')ijj{k2, fc)}*+(A?i)*+(fc2)|mc > . (9) 
and thus its squared norm: 

F{k,k')(x J dkidk2\tlj{ki,k2)^p{k,k') 

+ ij{h,k)ij{k2,k') + ^P{ki,k')^{k2, k)\\lQ) 




FIG. 2: State of four atoms, scattered pairwise back-to-back. 
Atoms detected by detectors A and B in an arbitrary pair of 
directions have partners recoiling in the opposite directions 
within a certain width imposed by the uncertainty on total 
and relative momentum of the atoms. The quantum state of 
the detected pair is obtained by a partial trace over the re- 
coiling momentum components, and the coincidence counting 
yield in the detectors is just the product of the single detector 
count signals 



This is the main result of the paper. Dealing explic- 
itly with the four atom component it is easy to see what 
happens. In Fig. 2. we illustrate the case of detection of 
a particle pair in random directions. Because, as noted 
below Eq. (4), ip{k, k') vanishes unless k and k' are anti- 
parallel, the first term in Eq. (10) only contributes if the 
detectors correspond to opposite directions. For oppo- 
site or random directions such as in Fig. 2, there is also 
no cross term between the second two terms because one 
vanishes whenever the other is finite. If {k, k') are nearly 
parallel, as illustrated in Fig. 3, the last two terms eval- 
uate the two different i/i-terms at the detector directions 
and at the direction specified by the integration variables. 
This means that values of the integration variables ki , k2 
exist (opposite to the detector directions), where both of 
the last terms in (|10[) contribute, namely if ki and k2 are 
within the "recoil cone" of both detection directions k 
and k'. For identical k and k' this interference give pre- 
cisely the factor 2 increase of coincidences compared to 
the case of random directions. We also note that the en- 
hanced coincidences occur within a solid angle specified 
precisely by this " recoil cone" . The next section develops 
our model one step further and carries out the calcula- 
tion for the special choice of a Gaussian Ansatz for the 
function ^(A:i,/c2). 
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FIG. 3: State of four atoms, with two atoms detected in nearly 
parallel directions. The detected atoms are not unambigu- 
ously identified with their recoiling partners if their momen- 
tum distributions are wide enough to overlap. This leads to 
an interference term in the coincidence counting yield in de- 
tectors A and B as expressed by the formal expression (I10|l . 



V. RESULT FOR A SIMPLE ANSATZ FOR 

ip{ki,k2) 

As stated above, the Bogoliubov approximation to the 
actual state can be found numerically by solution of lin- 
ear wave equations. In this section, we shall rather take 
a simpler approach by making an Ansatz for the shape 
of the function ip by appealing to the dynamics and the 
conservation laws valid in the bipartite collision dynam- 
ics. 

Energy conservation, which is effectively enforced dur- 
ing the temporal solution of the Schrodinger equation, 
suggests that the atomic pair state ip{ki,k2) describes 
particles with the same energy as the incident condensate 
particles. Momentum conservation suggests that they 
also have the same total momentum as the colliding pair. 
Due to the finite size of the colliding clouds, this does not 
strictly imply that the two atoms must have exactly op- 
posite momenta. The finite size of the collision zone im- 
plies a quantum mechanical momentum uncertainty, and 
if contributions from small regions (with correspondingly 
large total momentum uncertainty) are added coherently 
the resulting phase matching condition is not sharp if the 
entire collision zone has finite size. There are therefore 
quantum fluctuations of both the modulus and direction 
of the momenta. When the particles escape from the col- 
lision zone as illustrated in Fig. 1., they are also repelled 
by the mean field interaction with the two condensate 



components. Here we will not try to describe accurately 
the effect of this interaction, which is in any case small 
when atoms leave a condensate after receiving an initial 
kinetic energy large compared to the chemical potential 
(2^ . Note, however, that such mean field repulsion has in 
fact turned out to be tractable in the atom laser output 
from a condensate, where a generalized ABCD matrix 
formalism yields and analytical description of the propa- 
gation di. 

Define coordinates such that the nearly parallel k, k' 
of interest are close to the negative z-direction. Their 
partners at fci , k2 must both be close to the positive z- 
direction, and we shall assume the z-coordinates to be 
equal and opposite and only look at their x and y com- 
ponents. Their widths are related to the wave functions 
of the colliding condensates, both due to the amplitude 
of collisions out of these condensates and due to the mean 
field repulsion, and they are thus in general anisotropic. 
We restrict for simplicity the integration to one trans- 
verse coordinate (putting vector arrows on the arguments 
will yield the 2D result), and we assume that a single pair 
is described by a wave function, where the wave function 
amplitude for the recoiling partner has a bell shaped pro- 
file, that we for simplicity approximate by a Gaussian, 
centered at minus the coordinates of the detected parti- 
cle. The width of this wave function is parametrized by a 
momentum width K which thus represents both the mo- 
mentum width of the colliding condensate particles and 
the acceleration due to the mean field. 

Assuming thus the last two terms in (fTO|) to be of such 
Gaussian shape, and ignoring the first term which van- 
ishes for the geometry studied, we can explicitly calculate 
the coincidence signal: 



F{k, k') cx 




I exp(-((fci + fc)2 + {k2 + k'f)l2K^) 
+ exp(-((A:i+fc')' + (fc2 + fc)')/2if')|' 
oc l-hexp(-(fc-fc')V2i^^), (11) 

where we recall that k and k' here refer to (small) trans- 
verse coordinates of the detector directions with respect 
to a given axis, i.e. (fc — k') is the radial momentum 
of the outgoing particles multiplied with their mutual 
(small) angle in radians. 

We recover the Hanbury Brown and Twiss correlations, 
and we observe that the correlations persist for final state 
momenta within a distance from each other of the or- 
der of the quantum mechanical uncertainty of the total 
momentum of the atom pair escaping the collision zone. 
This is in accord with our interpretation in terms of the 
interference between the indistinguishable components il- 
lustrated as the overlapping recoil cones in Fig.l c), that 
leads to the last term in pT|) depending on both k and 
/c', whereas the direct terms lose the k and k' dependence 
due to the Gaussian integrals. 

It is interesting to note, that Ea. ([TT|) follows from 
a two-state amplitude il!{k,k') oc exp(— (A: -I- k'Y /2K^) 
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for transverse momentum components of atoms propa- 
gating in nearly opposite directions, and therefore the 
two atom component of Eq.® predicts a correlation 
of atoms in opposite directions with the dependence 
|V'(fc,fc')P oc exp(-(fc + k'fjK'^). The Hanbmy Brown 
and Twiss bunching thus occurs within a Gaussian width 
that is \/2 times larger than the range of correlation of 
recoiling atomic momenta. This prediction for the Gaus- 
sian wave functions has been verified by more detailed 
analysis of the full 3D propagation, [5| . One way to un- 
derstand the broadening is to recognize that the density 
dependence of the pair production mechanism results in 
a source which is spatially narrower than the condensates 
themselves. 

Although we have based our analysis on the four-atom 
component of the full many body states, we have argued 
that a calculation based on the full state would yield the 
same results, and in particular that the HBT correlation 
amounts to a factor of two in parallel directions while the 
correlation in opposite directions is not limited by this 
factor. When multiple scattering is neglected, the com- 
plete many-body problem is solved by the Bogoliubov-de- 
Gennes equations, and the resulting Gaussian/thermal 
character of the many-body state is fully accounted for by 
the second moments. This does not imply, however, that 
one would get the same quantitative results for scattering 
of few and many atoms. If the normalized wave function 
for a single scattered pair in the case of a low scattering 
probability is denoted x(fci, fc2), it may be reasonable to 
describe the collision process by the effective Hamiltonian 
H = K J dkidk2x{ki, k2)'^^ {ki)'^^k2) + h.c, where k is a 
coupling strength. The unitary time evolution operator 
is the exponential of this operator multiplied by (t/ih) or 
integrated over a suitable time interval. We note that this 
does not generally result in an expression for 'ip{ki, ^2) in 
(HI which is proportional to x{^l^ ^2)- In the case of sin- 
gle mode squeezing, it is known that one must evaluate 
the hyperbolic tangent function of the squeezing parame- 
ter to convert the squeezing operator to the normal order 
form and in our general multimode case, normal or- 
dering is accomplished by evaluating the tanh function 
of a matrix argument [l2l . [isj . For small arguments, in 
the perturbative regime of spontaneous four wave mixing, 
tanh is a linear function, and we get the same momen- 
tum dependence. Outside the perturbative regime, we 
retain the factor 2 bunching effect by our general argu- 
ment, but the precise shape of the correlation peak may 
be modified. 



VI. DISCUSSION 

We have presented a simple interpretation of the ob- 
served Hanbury Brown and Twiss correlations observed 
in the elastic scattering of Bose-Einstein condensates. We 
emphasize, that in order to make quantitative predic- 
tions, it is necessary to make a more elaborate calculation 
of the pair formation and the propagation of the atoms 



both in free space and in the regions where the mean 
field of the condensate components act as a perturbing 
potential. Such a description is offered by the Bogoli- 
ubov theory in Refs.[i, i, [l^, and we note that 10] as 
well as [4] also provide numerical evidence for the density 
correlations discussed in the present paper. Our interpre- 
tation relies on the structural property of the solutions to 
the Bogoliubov theory ([?]), but it proceeds by applying 
a different physical reasoning which recognizes that the 
two detected particles are accompanied by collision part- 
ners propagating in the opposite directions, and we hence 
observe part of a four-atom state. This is an appealing 
picture, in particular because the prediction of the coin- 
cidence signal, and in particular its width, relates to the 
transverse spreading of the pair wave functions of oppo- 
sitely propagating atoms after the bipartite collisions. 

As we discussed in the text, when observed from only 
one direction, the reduced density matrix of the expand- 
ing atomic cloud is similar to a thermal state. This den- 
sity matrix is sufficient to predict the outcome of any 
measurement on the observed part of the system, and it 
explains the experimental findings as an analogue of the 
observed bunching of the photons from a thermal/chaotic 
light source. The optical Hanbury Brown and Twiss ex- 
periment has a characteristic transverse spatial scale over 
which the correlation falls to unity, related to the trans- 
verse momentum distribution of the photons impinging 
on the detector, and in a similar manner we have a finite 
transverse coherence length in the atomic scattering ex- 
periment. We discussed the isotropic s-wave scattering, 
with possible corrections due to anisotropy of the collid- 
ing clouds and spatial phase matching. In addition, one 
may apply a periodic background potential, which may 
alter the energy dependence on the momentum vector of 
moving atoms, and hence modify the scattering profile 
[2^ |25||. and with confinement to one dimension, it may 
lead to highly selective population of specific momen- 
tum states with strong, observable quantum correlations 
[2i,[23. 

It is interesting to recall that a mixed quantum state, 
i.e., a density matrix, can always be formally obtained as 
the reduced state of a larger quantum system which is in 
a pure state, and in particular any thermal quantum state 
of a bosonic degree of freedom can be modelled by a pure 
squeezed state in a doubled tensor space. This is known 
as the "thermofields" formulation [2g|, and for example 
the single mode thermal state ([2|) can be obtained as the 
trace over one of the modes of a non-degenerate two mode 
squeezed state, as obtained, e.g., from a non-degenerate 
optical parametric oscillator (OPO), 

IV-OPo) (12) 

n 

In our four- atom analysis, the apparently thermal state 
arriving at nearby detectors, is precisely part of such a 
larger system. The advantage of this insight is that the 
spatial scale of the extended state, in our case the prob- 
ability distribution of the total momentum of scattered 
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atoms, directly yields the density correlations in the re- 
duced density matrix. 

Finally, if the collision occurs between two different 
bosonic species, the same kind of correlations will oc- 
cur for the density correlations of each species, but not 
for the cross correlation, where the recoiling atoms are 
distinguishable, and hence do not interfere. For colli- 
sions between bosons and fermions the situation is dif- 
ferent. Electrons have been demonstrated to show anti- 
bunching related to their fermionic character (2^ . and 
anti-bunching has also been demonstrated for neutral 
fermionic atoms [s^, [3l| . Collisions between a Bose con- 
densate and a degenerate Fermi gas, where all fermions 
initially occupy orthogonal states, but where Pauli block- 
ing forbids more than one atom ending up in the same 



final state should lead to observable anti-bunching effects 
in the scattered bosons. 
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